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ABSTRACT 

The physical parameters of the stars in the central region of the globular cluster 
M5 (NGC 5904) were determined from UBV photometry using Kurucz (1992) (K92) 
synthetic flux distributions and some empirical relations. It is found that the bluest 
horizontal branch (HB) stars have higher luminosities than predicted by canonical 
zero-age horizontal branch (ZAHB) models. Parameters of the mass distribution on 
the HB stars are determined. It is shown that the gap in the blue HB previously 
reported (Markov et al. 1999, hereafter Paper I) is probably a statistical fluctuation. 

Key words: globular clusters: general - globular clusters: individual (M5=NGC5904) 
- stars: evolution - stars: horizontal branch - stars: Population II. 
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1 INTRODUCTION 

Crocker, Rood & O'Connel (1988, CR088) determined spec- 
troscopically logTcff and logg for 14 BHB stars in globular 
duster (GC) M5. They found that M5 BHB stars matched 
the canonical ZAHB models quite well. In our photometry 
(Paper I) there are 11 stars common with CROSS. Assuming 
a relationship between the temperatures and surface grav- 
ities as derived in CR088 and the atmospheric models of 
K92, we are able to obtain the physical parameters for more 
than 100 BHB stars from our photometric observations. 

According to Buonanno et al. (1981) (BCF81) and 
Crocker & Rood (19S5) the bluest BHB stars tend to rise 
above the canonical ZAHB whilst Bohlin et al. (1985) did 
not find such effect. The completeness of our sample throws 
additional light on this matter. 

The present work is structured as follows. In the 
transformation of observed colours and magnitudes to the 
(logTeft, logL) plane is described. In an approximation 
of the optical colour-magnitude diagram (CMD) in terms of 
synthetic HBs is presented and the statistical importance of 
the gap in the blue HB discussed in Paper I is evaluated. §^ 
is dedicated on the position of BHB stars in the theoretical 
diagram and finally our findings are summarized in §H 



2 PHYSICAL CHARACTERISTICS OF HB 
STARS 

We determine the effective temperatures of BHB stars from 
their reddening-free colours. For this purpose, we use theo- 
retical models of energy distributions for different tempera- 
tures, gravities and abundances, in particular the line blan- 
ket atmospheres of K92. We choose the model with metal- 
licity [Fe/H]= —1.5, i.e. the nearest value to that derived in 
Paper I ( — 1.38). It is well known that the model grids make 
folds on the {U - B)o, {B - V)q and {U - V)o, {B - V)o 
planes. This leads to a nonsingle- valued transformation from 
colours to temperatures and gravities. Unfortunately, all the 
BHB stars of M5 are situated just within the folds and we 
forfeit the opportunity for simultaneous determination of 
logTcff and logg. After a suitable gravity-averaging of the 
atmospheric parameters, it becomes possible to achieve the 
necessary unequivocal dependence of the colours on the ef- 
fective temperature. 

The analysis of the data for BHB stars is divided into 
two parts. We consider separately hot BHB stars (those with 
(B - V)o < -0.02) and cool BHB stars (with {B - V)o > 
—0.02). The reason is that the different colour indices as 
temperature and gravity indicators exchange their role at 
about Tofl = 9 200K. The {U — B)o colours as a measure of 
the Balmer discontinuity and {U — V)o as a long baseline 
colours are very sensitive temperature parameters for B- 
type stars in the UBV system. The {B — V)o colour is a 
good temperature indicator for cooler stars (log Teg < 3.96). 
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Figure 1. Comparison between the effective temperatures 
log Tpff derived in this work and the CROSS values for 11 common 
stars. 



For that reason different models are adopted to analyse the 
hot and cool BHB stars. 



2.1 Determination of logToji for hot BHB stars 

The effective temperature is derived from the observed UBV 
data using the relations Q(logroff) and {U — V^)o(log Tcff ) 
found in K92. These relations are applied for Q < 0.0 and 
{U — V)o < 0.10. Both relations are employed for logy in 
the interval (3, 4.5) and then the derived T^ff values are 
averaged. The effective temperatures for 67 BHB stars are 
determined in this way. The dispersion in log Tctt in the two 
cases is ±0.013 and ±0.020. The mean standard errors of Q 
and {U — V)q are about ±0.04 mag, so the error in log Tcff 
should not exceed ±0.043 (±1 lOOK at T^s = lOOOOK). 



2.2 Determination of logTcft for cool BHB stars 

Q and {U — V)o are no longer sensitive temperature indi- 
cators for stars in this evolutionary stage. {B — V)o with 
a maximum gravity dependence of A{B — V)o — 0.18 at 
log Tcff — 3.9 for logg in the interval (2, 4.5) is a more ap- 
propriate indicator. An assumption about logy is needed 
to derive the temperatures. Model computations show that 
{U — B)q is a strong gravity indicator for these stars. This 
is demonstrated in Fig. 7 of Paper I, where the theoretically 
predicted two-colour relations taken from K92 are overlaid 
on the observed data. As we discussed in §4 of Paper I, the 
majority of these IfB stars show an ultraviolet deficiency 
and their positions are consistent with the loci of the mod- 
els with log (J ~ 2.0. For this reason logg within the interval 
(2, 3.5) is chosen here. The model suggests a mean standard 
deviation for log Tcff of ±0.021. Taking into account that 
the uncertainty in {B — V)o is about ±0.04 mag, we find 
the error in logTcfi ~ ±0.045, which corresponds to 900K 
at Teff =8 OOOK. 

The effective temperatures of 111 hot and cool HB stars 
are obtained using this technique. A check of the validity of 
the determined temperatures is made by comparison with 11 
BHB stars in common between the CROSS and our obser- 
vations. Figure |l] illustrates a tendency of our temperatures 
to be higher on average by 0.008 ± 0.003. 

Provided that the effective temperatures are de- 
termined, appropriate K92 relations (BC(logrcfi) and 
BC [{B — V)q] for hot and cool BHB stars, respectively) 
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Figure 2. (logL, logTefi) diagram for all 111 BHB stars (solid 
circles) with UBV magnitudes in our sample. The open triangles 
denote 14 stars of CROSS. The solid hne is the ZAHB of Dorman 
et al. (1993). 



are used to derive bolometric corrections (BC), and then 
bolometric magnitudes and luminosities. The derived lumi- 
nosities are based on the assumption that (m — M)vq ~ 
14.46 ± 0.03 and M®; = 4.72. The best fit of our ob- 
served CMD to the canonical ZAHB of Dorman et al. 
(1993) (DR093) is used to obtain the distance modulus. 
The DR093 evolutionary models are for [Fe/H]= -1.48, 
[0/Fe]= 0.63, Yhb = 0.25 and core mass A'/corc = 0.485 Mq 
(hereafter we imply only this set of DR093 tracks when 
we refer to the evolutionary theory). The mean square er- 
ror of the luminosity determinations is ±0.1 in logL. The 
(log L, log Tcff) diagram for all 111 BHB stars with UBV 
magnitudes in our sample is given in Fig. ^. The stars of 
CR088 with spectroscopically measured atmospheric char- 
acteristics are also plotted as open triangles. 



2.3 Determination of the stellar masses 

Knowing the effective temperatures, luminosities and sur- 
face gravities of the stars, we can estimate the corresponding 
masses independently of their evolutionary history through 
the relation M/L ~ g/T^g or in more convenient form 

log (M/Mq) = log (L/Lq) - 4 log Tcff ± logy ± 10.611. 

We rely on a sample of 12 BHB stars with reliable log g 
and log Toff data (CR088) to derive a linear relation between 
surface gravity and effective temperature. A least square fit 
yields: 



logy = -11.244(±1.929) ± 3.699(±0.478) logTrf 



(1) 



Equation (|l]) is deduced for stars with log T^g in the interval 
(3.948, 4.145), and it is used to obtain the surface gravities 
of all BHB stars with r.m.s. ~ 0.15, which corresponds to 
the errors in the input data. 

The errors of the mass determination consist of both 
the errors of log T^s and log y , and those of log L. Assuming 
errors in logTes of ±0.043, in logy of ±0.15, and in logL of 
±0.1 we obtain a mean error in the mass derivation of ap- 
proximately ±0.19. Changes in the adopted distance modu- 
lus and extinction would result in a systematic shift of the 
derived masses. The dependence of the stellar mass on the 
effective temperature is shown in Fig. ^ The distribution 
of the masses of 62 stars situated within the temperature 
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Figure 3. (mass, log T^ff ) diagram. The solid circles are stars with 
obtained masses. The 12 stars of CROSS are plotted as open trian- 
gles. The solid line is the (ZA mass, logTejf) relation; the dotted 
line is the relation between the stellar mass and maximal tem- 
perature reachable during the post-ZAHB evolution; the dashed 
line is the isochrone at t = 110 Myr. The inset shows the mass 
distribution. 



range (8 800K, 14 000 K) of validity of the linear relation 
logg (logTcff ) is plotted in an inset on Fig. ^. 

During their HB lifetime, the stars reach a maximal 
effective temperature (see, for instance, DR093). This tem- 
perature as a function of stellar mass is also represented in 
Fig. ^by dotted line. It almost completely coincides with the 
isochrone at t — 60 Myr. Although the error of the stellar 
masses is considerable. Figure ^ implies that the masses of 
the hottest stars (logTeff > 4.07) are systematically higher 
than those corresponding to the limit temperatures. This 
suggests a discrepancy in luminosity compared with theory. 
The masses of the 12 stars of CR088 are also plotted in Fig- 
ure ^ (open triangles) . The big scatter of these stars around 
the zero-age relation can be interpreted as an inaccurate de- 
termination of the surface gravities (see Table 2 of CROSS). 
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Figure 4. Colour-magnitude diagram of the HB region of the 
M5 and the ridge line (dashed line) used to build the stellar dis- 
tribution along the HB. Solid line is the ZAHB of DR093. 
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3 HB SYNTHETIC (My, {U - V)o) DIAGRAM 
3.1 Parameterisation of the HB 

A curve (ridge line), which is close enough to a HB strip, is 
determined by means of averaging the stellar positions on 
small segments covering the HB. A position along the HB, 
from the red extreme toward the blue end, is given by arc 
length (^hb) of the normal projection on the ridge line (for 
similar procedures, see Rood and Crocker 1989; Ferraro et 
al. 1992; Dixon et al. 1996). 

The ridge line and HB region of M5 is shown in Fig. ^ 
(dashed line). To make possible a comparison of our distri- 
bution to other ones, we give the scale ratios A^hb/A(C/ — 
V)o = 3.341 and AIkb/AMv = 2.468. The values of the 
HB red extreme of the ridge line are My^ = 0.43 and 
{U — V)^^ = 0.7, the blue and red edges of the instabil- 
ity strip at the level of the ridge line are {U — VY^^ — 0.26 
and {U - V)f^ = 0.45 (see Fig. §). The observed distribu- 
tion is shown in Fig. ^ both as a histogram and a cumulative 
frequency. The ^HB-values of the blue and red edges of the 
instability strip are £^^^ = 0.244 and = 0.143. 



Figure 5. The observed stellar distribution along the HB in M5. 
The histogram of the £HB~coordinate of the stars is represented 
on panel (a) and the cumulative frequency - on panel (b). 



3.2 Theoretical models of the HB 

To create a model HB, we transform the mass-age rela- 
tion (after ZAHB) into a colour-magnitude plane using the 
theoretical correspondence between the (mass, age) and 
(logL, log Toff) or (logg, logTcff) relations. Evolutionary 
tracks (DR093), synthetic atmospheres (K92), the adopted 
metallicity, and distance modulus are the same as in A 
Hermite interpolation (Kahaner et al. 19SS) of the DR093 
tracks is employed. We adopt both a unimodal (Gaussian) 
mass distribution and a bimodal one (see e.g. Rood 1973; 
Lee et al. 1990; Dixon et al. 1996). 

An assembly of 76 RRLyrae stars randomly distributed 
in the instability strip in (My, {B — V)o) plane with a mean 
magnitude (Vrr) — 15.09 (Reid 1996) is added. From the 
CMD, the blue edge of the instability strip is located at 
{B - V)o^^ = 0.16 and the red edge - at (B - V)f'' = 
0.38 (for comparison see BCFSl). The adopted number of 
variable stars is derived under condition that the sample 
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Moacl 


Parameter 


Permissible interval 


Best fit 


UMD 


(A/hb) 


(0.5801, 0.6782) Mq 


O.6327M0 




°'AfHB 


(0.0179,0.0433) Mq 


O.O321M0 


BMD 


(Mhb>i 


(0.5699, 0.6432) Mq 


O.6O57M0 




'^(Mhb),1 


(0.0181, 0.0199) Mq 


O.O192M0 




(■WhB>2 


(0.6110, 0.6850) Mq 


O.6493M0 




'^(A/hb),2 


(0.0258, 0.0288) Mq 


O.O274M0 




/l 


(0.3703, 0.4054) 


0.3898 



Table 1. Permissible intervals of the parameters determining the 
model mass distributions and their best-fitting values. 

represents 74 per cent of all HB stars located in the measured 
region. 

Taking into account the weak dependence of the HB 
lifetime on the total mass of the stars, we can introduce 
a universal HB lifetime 110 Myr. Due to the assump- 
tion of a constant rate of populating the ZAHB, the ages 
after the ZAHB are considered as random numbers uni- 
formly distributed in the range 0-110 Myr. We also presume 
known stellar mass distribution as well as age distribution. 
For a unimodal mass distribution (UMD) the transforma- 
tion method (Box-MuUer algorithm) is used to produce the 
Gaussian deviates, whereas in the case of a bimodal mass 
distribution (BMD) the appropriate deviates are generated 
by means of the rejection method. Normally distributed ran- 
dom 'observational' errors are introduced at the end of the 
mapping. 

To restrict the unknown parameters of the mass distri- 
bution that produce the observed CMD, we compare statis- 
tically the corresponding distributions of synthetic and ob- 
served HBs. A two-sample Kolmogorov-Smirnov (KS) test 
is used to check the null hypothesis that model and ob- 
served cumulative ^HB-frequencies (^HB-curves) are drawn 
from the same parent distribution. Simultaneously a chi- 
square (x^) test is applied to binned ^HB-coordinate distri- 
butions. Each synthetic HB contains the same number of 
stars as our observed sample. At the 99 per cent level of 
significance, we disprove the null hypothesis for all parame- 
ters outside the regions listed in Table I (Mhb) and 
are the mean value and the standard deviation of the Gaus- 
sian mass distributions, respectively. In the case of bimodal 
mass distribution, we have an additional parameter /i, the 
fraction of HB stars contributed by the first Gaussian. The 
mean values of the significant levels (Qks) and (Q^^) ^"-^ 
obtained by generating ~ 250 000 synthetic HBs with the 
best-fitting parameters (see Table hi). When the UMD is 
used (Qks) = 0.43 and {Q~^2^ — 0.17, whereas in the case 
of BMD (Qks) = 0.48 and (Q^2) = 0.25, which makes the 
second distribution preferable. 

Figure ^ shows four synthetic HBs computed by means 
of UMD, as well as the corresponding cumulative distribu- 
tions of coordinates, the histograms of masses and temper- 
atures. Similar results are plotted on Figure ^ for bimodal 
mass distribution. In both cases (UMD and BMD) the cor- 
responding four synthetic HBs are chosen randomly from 
the whole sets of HB simulations. At first glance, there are 
considerable differences between the separate synthetic HBs 
although the model assumptions are just the same. The prin- 
cipal cause is that the random deviates used for the mass 




Figure 6. Four synthetic HBs randomly chosen from the set 
of ~ 250 000 HB simulations and cumulative frequency of £uB^ 
coordinate, which are computed by unimodal mass distribution. 
The insets show the histograms of the T^ff and mass distributions. 
The dashed line is the ridge line. 

generation only asymptotically approach the desired distri- 
butions. However, a thorough inspection reveals that the 
blue tails of UMD synthetic HBs are less populated and 
more extended than tails of the BMD synthetic HBs - a 
prospective effect due to the bigger distribution tail area of 
the unimodal mass distribution. 

The main discrepancy between the observed and model 
HBs is the downward deviation of the synthetic blue tails 
from the ridge line. Both UMD and BMD computations 
show the same tendency. The formal cause may be the incon- 
sistency of the uniform distribution of the stellar ages after 
ZAHB adopted in the computations, although, most proba- 
bly, the real cause is the unreliability of canonical evolution 
theory. 

3.3 Gaps on the blue HB 

Even a passing glance at the simulated HBs ascertains sub- 
stantial statistical fiuctuations, which modify their morphol- 
ogy. The observed gap discussed in Paper I may thus be 
such a fiuctuation. We consider as gaps regions with a local 
minimum of the slope of a monotonically increasing cumula- 
tive frequency function. The standard procedure to test the 
statistical significance of suspected gaps is based on the as- 
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Figure 7. Same as in Figure ^, but for the bimodal mass distri- 
bution. 



sumption that the probabihty to find stars along the cumu- 
lative distribution is uniform with a value Pl left of the gap, 
and uniform with another value Pr right of the gap (Aisen- 
man et al. 1969). Usually, the significance of this hypothe- 
sis is not taken into account. Instead, one verifies whether 
the probability of finding stars in the gap substantially dif- 
fers from the mean of the probabilities Pl and Pr (Howar- 
den 1971). When applied to cumulative ^HB-distributions of 
HB stars, the standard method yields a local highly over- 
estimated statistical significance of the gaps. In the present 
work, we apply KS-test to a section of i^HB-curve in length 
of five gap widths in order to examine the hypothesis of 
uniform distribution around and in the gap (probably with 
three different densities). 

The observed gap should be represented as a domain 
on the {W, Ins) plane {W is the gap width) because of 
photometric errors. To derive the borders of this region, we 
vary the stellar magnitudes and colours within photomet- 
ric accuracy bounds reported in Paper I, and produce a set 
of ~ 200 000 modified observed HBs. Then, we obtain the 
joint probability density function (PDF) of the observed gap 
fo {W, Ins) using an automatic search and inspection of aris- 
ing intervals. Three elliptic-like contours (with major axes 
in ordinate direction) containing 95, 50, and 5 per cent of all 
detected gaps are plotted on Figure pl to depict this function. 
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As mentioned in 5)3.2, the difference between the mod- 



Figure 8. Characteristics of the observed and detected gaps on 
the synthetic HBs computed for UMD and BMD. 



els and observed data is statistically non-significant. This 
enables us to use the estimation of a mere chance proba- 
bility that a gap comparable with the observed one may 
appear on the synthetic HBs. We consider the width and 
the location of the biggest gap on the simulated £HB-curves 
as random variables. The significant level Qks of the gap re- 
gion can also be interpreted as a random parameter, which is 
a function of the gap sharpness. We construct the sampling 
distributions empirically by using the set of ~ 250 000 sim- 
ulated HBs. The joint PDFs fs (WJhb) and (QksJub) 
are depicted on Figure |^ as isodensity contours containing 
95, 50, and 5 per cent of all revealed gaps. The gaps are also 
plotted as points. 

To test the null hypothesis that there is no genuine dif- 
ference (within the limits of random variability) between the 
parameters of observed and model gaps, we use the follow- 
ing three probabilities (see Appendix Pi and P2 are the 
averaged (over fo) probabilities of arising of synthetic gaps, 
which are wider than the observed one. They are calculated 
by the marginal PDF {ms{W) for the (W, £hb) relation 
and ms{QKs) for the {Qks, Ihb) relation) and by a condi- 
tional PDF (/s(W^|/hb) and fs{QKs\£iiB)), correspondingly. 
P2 takes account of location on the ^HB-curve as distinct 
from Pl. P3 is the probability of appearance of synthetic 
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UMD model 


BMD model 


Relation 


Pi 






Pi 


P2 P3 


{W, Ihb) 


0.21 


0.10 


0.30 


0.20 


0.20 0.48 


(log Qks-, (-hb) 


0.08 


0.08 


0.16 


0.09 


0.08 0.17 



Table 2. Probabilities used to test statistical significance of the 
observed gap. 

gaps outside the region enclosed in the averaged (over /o) 
isodensity curve on which the observed gap is situated. The 
so-defined probabihties for the (VK, ^hb) and [Qks, ^hb) 
relations are given in Table |^. By the conventional standard 
of statistical significance (probability < 0.05), the observed 
gap is non-significant. Therefore, the observed gap is not so 
extraordinary, as it seems to be, and we tend to consider it 
as a statistical fluctuation. 



4 BHB STARS AND THE THEORETICAL H-R 
DIAGRAM 

In the (logL, log Toft) plane (Fig. |^, the tendency for the 
blue end of the IfB to rise above the ZAHB (the oxygen- 
enhanced models (DR093) for [Fe/H]= —1.48) is obvious. 
BCF81 and Crocker & Rood (1985) have found the same 
effect, using different methods and models for transforming 
the observational data from a (V, {B ~ V)) plane into the 
theoretical (log L, log T^s) plane. On the other hand, Bohlin 
et al. (1985) have not confirmed this trend (their Figure 7) 
when perform ultraviolet observations with a rocket-borne 
telescope near 1540A and 2360A for 50 BHB stars. 

An important argument for the validity of an upward 
trend is the location in the (log L, log T^s) diagram of those 
14 BHB stars with atmospheric parameters being already 
measured spectroscopically (CR088). They are plotted in 
Fig. ^ as open triangles. On canonical theory we expect the 
observed HB to be close to the ZAHB (DR093). One pos- 
sibility for the revealed disagreement is that we have un- 
derestimations of the adopted photometric errors. To check 
this we use 11 stars in common with CR088. From the 
K92 temperature-colours calibrations we compute the cor- 
responding {U — V)o, {U — B)q and {B — V)o values for 
each star. The mean differences between the model and 
the observed colours are as follows: A{U — V)o ~ 0.02, 
A{U - B)o = 0.01, and A{B - V)o = 0.02 - values compa- 
rable within la of our photometric errors. 

It can be seen from Fig. ^ where the theoretical ZAHB 
is fitted to our observed data, that the blue end of the ob- 
served HB (dots) do not match the model of DR093 (solid 
line). This discrepancy has already been noted by Sandquist 
et al. (1996). The fiducial line of BHB derived by these au- 
thors is plotted for comparison as dashed line in Fig. |^. 
Sandquist et al. (1996) assume that the disagreement re- 
sults from the colour calibration. However, the fact, that 
the high temperature stars with spectroscopically obtained 
temperatures also have much higher luminosities (see Fig. ^ 
than the predicted ones, means that the observed effect can- 
not be explained by the colour calibration only. In Fig. ^b, 
we compare the positions of these stars (open triangles) 
(the V mag for three of them are taken from BCF81) with 
the location of the model HB. It is evident that stars with 
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Figure 9. (My, {B-V)o) diagram (panel (a)) and (My, logTcff) 
diagram (panel (b)) for BHB stars with UBV photometry. Open 
triangles denote the stars of CR088. Solid line is the ZAHB 
(DR093). Dashed line on the panel (a) is the fiducial line of 
Sandquist et al. (1996). 



logTcft > 4.04 deviate upward from the model with AMy 
in the range (—0.2 to —0.5). These values are too large to be 
explained with photometric errors. A similar effect of 'over- 
luminosity', reaching ~ 0.4 mag has been found by Grundahl 
et al. (1998) in the globular cluster M13 using Stromgren 
photometry. Moreover, according to Grundahl et al. (1999) 
all stars with Tch > 11 500 K in 14 GCs demonstrate a sys- 
tematical deviation from the canonical ZAHB in the sense of 
being brighter or hotter than the models. This effect is inde- 
pendent of the metallicity, the mixing history on the RGB, 
and other cluster parameters. The results in Fig. ^ are in a 
remarkable agreement with their results. 

In explanation of this effect, we note that apart from 
the most popular picture, in which HB stars are distributed 
along ZAHB because of differences in the amount of mass 
loss during the 'He-flash' at the RGB tip, Fusi Pecci & Ren- 
zini (1978) focused attention on the possible effect of stellar 
rotation. With a suitable choice of parameters the blue end 
of the 'rotational' ZAHB is expected to be more luminous 
than the non-rotating canonical one. 

The standard method used in §^ to obtain physical pa- 
rameters of the BHB stars does not permit simultaneous 
determination of the effective temperature and the surface 
gravity. The latter parameter is lost during the inevitable 
averaging inherent to the elimination of the ambiguity of 
the Kurucz grids (see §^). For that reason we cannot con- 
firm independently the effect found in 6 GCs including M5 
(Moehler et al., 1995; Moehler, 1999), which consists in a 
displacement from the ZAHB to lower logg of the bluest 
HB stars having spectroscopically measured surface gravi- 
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ties. We shall only remark that the linear relation (|l|) shows 
the same tendency to 'under-gravity' towards the bluest HB 
stars in M5. Grundahl et al. (1999) reveal the same connec- 
tion between the 'over-luminosity' and the 'under-gravity' 
in 14 GCs, and suppose that these effects are different man- 
ifestations of a uniform physical mechanism. Relying on the 
spectroscopic data for field and GC BHB stars, the authors 
(Grundahl et al. 1999) hypothesize that this physical mecha- 
nism consists in radiative levitation of elements heavier than 
carbon and nitrogen into the stellar atmosphere, rather than 
a stellar interior or evolution effects. 

Formally, the possibility for a displacement toward 
higher luminosities of some stars due to evolution off the 
ZAHB cannot be ignored. Moreover, when the BHB stars 
reach the HB lifetime they demonstrate the 'under-gravity' 
effect simultaneously with the 'over-luminosity'. However, 
according to the canonical evolutionary theory (DR093) 
and the adopted concept of constant evolution rate onto the 
ZAHB, there is no well-grounded reason that so many low 
mass stars being at the same time in a transient evolutionary 
phase as core helium exhaustion. 



5 SUMMARY 

In the present paper we have applied standard methods to 
determine the stellar physical parameters using a new UBV 
stellar photometry for the HB stars of the GC M5. Com- 
paring the observed UBV colours with K92 synthetic at- 
mospheres, we have derived the effective temperatures and 
luminosities of a large sample of BHB stars. The surface 
gravities have been derived assuming a linear relation be- 
tween spectroscopically measured log g and Tcff of 12 stars 
(CR088). 

We have found considerably high luminosities for the 
bluest HB stars than predicted by the canonical theory. 

Approximations of the distribution of stars in the 
(My, [U — V)o) diagram have been performed under as- 
sumptions of unimodal and bimodal mass distributions. 

The statistical importance of the observed gap in the 
blue HB at ^hb = 0.42 has been studied by simulations. We 
conclude that it is a statistical fluctuation. 
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APPENDIX A: DEFINITIONS OF THE 
PROBABILITIES USED TO TEST 
STATISTICAL SIGNIFICANCE OF THE 
OBSERVED GAP 

For the sake of completeness, we present definitions of the 
probabilities Pi, P2, and P3 used in §3.5. We restrict our 
attention only to the (W, £hb) relation. 



P2 = 




(W) dW ruo (W) dW , 



1- J .fs{W,£HB)dW) 

J_ 

rus (fee) 



where /o (PV, i^Hs) and fs (W, ^hb) are joint PDFs of the ob- 
served and synthetic gaps, rris iW), rria {£hb) and irio (W) 
are the marginal PDFs of the fg (VF, ^hb) and fo {W,£b_b)- 
The integration is performed over domain in {W,£ub) plane 

with fo{W,£HB) > 0. 



P3 = y y fsiWJuB) dWd£HB 

with a region of integration determined by the inequality 
fs{W,£iiB)< I I fs {W, £hb )fo(W, £nB) dW dfes . 
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